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A full angular analysis has been performed for the pseudo-scalar to vector-vector decays JS" ^ 
J/tp K*^ and J/tp cj>to determine the amplitudes for decays with parity-even longitudinal (Ao) 

and transverse (Ay) polarization and parity-odd transverse {A±) polarization. The measurements 
are based on 190 B" candidates and 40 -Bg candidates collected from a data set corresponding to 
89pb~^ of pp collisions at -^s = 1.8 TeV at the Fermilab Tevatron. In both decays the decay 
amplitude for longitudinal polarization dominates : \Ao\'^ — 0.59 ± 0.06 ± 0.01 for the decay 
and l^oP = 0.61 ± 0.14 ± 0.02 for the B^ decay. The parity-odd amplitude is found to be small: 
l^xP = 0.13+° Qg ± 0.06 for B° and \A±\'^ = 0.23 ± 0.19 ± 0.04 for B°. 

PACS numbers; 13.20.He, 13.25.Hw 

The decays J /ip K*^ and J/tp (p are pseudo-scalar to vector- vector decays and in principle have three 

decay amplitudes which can be determined by studying the angular distributions of the final state particles. These 
decays can have orbital angular momenta between the J /ip and K* (or 0) of 0, 1, or 2, and three decay amplitudes 
govern these transitions. Another convenient description is given in the transversity basis in which the decay 
amplitudes are defined in terms of the linear polarization of the vector mesons. In this basis the L = 1 (P wave) 
decays are governed by a single decay amplitude, corresponding to a parity-odd correlation between transversely 
polarized vector mesons. The other two decay amplitudes, f-nd , are combinations of the parity-even L — Q and 
L = 2 (S and D wave) decays. The longitudinal polarization fraction, as is commonly defined in the hclicity basis 
is given by Tl/T = |Aop. 

Determination of the decay amplitudes and phases probes the limitations of the factorization hypothesis Fac- 
torization assumes that a weak decay matrix element can be described as the product of two independent (hadronic) 
currents, in this case treating the J /ip and B K* {(j>) as currents. To the extent that factorization holds, final state 
interactions are negligible. The matrix elements factorize into short and long distance (weak and strong) parts and 
the polarization decay amplitudes do not interfere and so are expected to be relatively real. 

A measurement of the parity-odd amplitude, A±, is pertinent to studies of CP invariance. For example, if the 
decay — s- J/ip K*'^ (with K*'^ — > tt'^) were to occur mainly through either a parity-odd or even amplitude, then 
this mode could be used as readily as the B^ J/%pK^ decay mode for determining |^ the CP nonconserving 
parameter sin(2/3). The situation holds as well for the decay B^ — + J/'^<f)i which is expected to have a very small 
CP decay rate asymmetry in the Standard (CKM) Model. In addition the polarization of the decay 5° J jip (p is 
interesting for its potential to improve the precision of measurements of the lifetime difference between the B^ mass 
eigenstates via an angular analysis 

The longitudinal polarization in — > J/^ K*'^ was first measured by ARGUS fflas ~1.0 followed by a CLEO ||] 
measurement of 0.8 ± 0.2 and a CDF (Collider Detector at Fermilab) measurement |1^ of 0.65 ± 0.11. A full angular 
analysis by CLEO ||ll| obtained a longitudinal polarization fraction of 0.52 ± 0.08 and a parity-odd fraction of 
0.16 ±0.09. The only previous measurement of the B^ decay mode, by CDF |l0), obtained a longitudinal polarization 
fraction of 0.56 ± 0.21, consistent with the results for 5° as expected under S'[/(3)-flavor symmetry. 

This Letter describes a full angular analysis of both the B^ and B^ decay modes based on 89pb~^ of pp collisions 
at ^/s =1.8 TeV collected with the CDF detector at the Fermilab Tevatron. In this analysis, the ^ J/ipK*^ and 
J/ipff) decays are reconstructed from the decay modes J/ip fi-~^fi'~, K*^ — > K'^ir^, and (p — > K~^K~. After 
the selection described below, there are 190 5° J/tpK*'^ and AQ B^ ^ J/ip(p candidates above background. The 
results presented here are independent of those in Rcf. . 

CDF is a general purpose detector and has been described in detail elsewhere [|l^. For this analysis important 
elements of the detector are the silicon vertex detector, with track impact parameter resolution of 10/ini, and the drift 
chamber, with charged particle momentum resolution of Sp^/pr ^ 0.001 p^i where px (in GeV/c) is the component 
of the momentum transverse to the pp collision axis. Chambers for muon detection provide coverage for particles with 
direction within 40 degrees of the transverse direction. 

The event sample for this measurement is selected by online trigger and offline reconstruction criteria. The trigger 
signature is the decay of a J/ip to two muons. Hits in the muon chambers forming a track segment consistent with 
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a nominal muon transverse momentum pT>3.0GeV/c satisiy the first level of the trigger. Candidate muon track 
segments must be separated by more than 10° in the plane transverse to the collision axis. The second level requires 
drift chamber tracks, with pT>2GeV/c, which extrapolate to the track segments in the muon chambers. The third 
level accepts J/'ip candidates with a reconstructed mass between 2.8 and 3.4GeV/c^. 

The offline reconstruction first selects J /ip candidates. Two muons of opposite charge satisfying quality requirements 
are combined into a J/?/; candidate. A kinematic fit requiring the muons to originate from a common vertex is 
performed, and the confidence level of the fit is required to be at least 0.01. This results in a sample of about 290,000 
candidate J/iJ^'s. 

Of these, candidates within 80MeV/c^ of the J/^ mass (3096.88 MeV/c^) ||l^ are then combined with two oppositely 
charged tracks that form a K* {(f>) for the B'^ — > J/tp K*^ {B^ ^ J/tp (f>) decays. CDF lacks the particle identification 
capability to distinguish between the two K-tt charge assignments for a K* candidate. However, the mass distribution 
for misidentified candidates is broader than the natural width of the K* and their contribution is largely suppressed 
by choosing the assignment yielding a K* mass closer to, and within 80MeV/c^ of, the world average (896.1 MeV/c^). 

To further improve the signal to noise ratio, the B'^ candidate is required to have pT>6.0GeV/c, with the K* 
having pT>2.0 GeV/c. All four particles from the B'^ decay are required to come from a common 'secondary' vertex, 
with the confidence level of the fit being greater than 0.001. The proper decay length (cx proper decay time) for the 
B'^ is required to be at least 100 //m. The resulting K-k^^ mass distribution is plotted in Figure |l|(a). 

For the B'^ decay, the (f) candidate is required to have a mass within lOMeV/c^ of the nominal (p mass 
(1019.4MeV/c^). The narrow natural width of the (p provides for better background rejection than the K* . In 
order to compensate for the lower expected yield of B^ events, looser requirements are imposed: transverse momen- 
tum cuts of 1.5 GeV/c and 4.5 GeV/c are applied to the (f) and B^ candidates. The confidence level of the four track 
fit is again required to be greater than 0.001, and the proper decay length is required to be at least 50 /im. The 
resulting KK^^ii~ mass distribution is plotted in Figure |^(b). 

The linear polarization of the vector mesons is determined from the angular distribution of the final state decay 
products. The decay angles are defined |^ as 0_r-*, 8t, and <i>T- In the rest frame of the K*, Qk* is the angle 
between the K momentum and the direction opposite the B meson (for the cj), the defines the angle). The 
transversity angles, Ot and $t, are defined in the J/ip rest frame: Ot is the angle between the /i"*" momentum and 
the perpendicular to the K*-K plane; <i>T is the azimuthal angle from the K* to the projection of the /i"*" momentum 
onto the K*-K plane. 

The decay angular distribution is ||^: 

dV/dn cx 2 cos2 Gk- (1 - sin^ Ot cos^ <^t)\Ao\^ 

+ sin2eK.(l-sin2eTsin2$T)|^j|P 

+ sin^e^. sin^exl^iP 

+ -i|Sin2eK. sin2eTsin2$TR.e(ASA||) 

=F sin^ Ok- sin26T sin$T Im (A| Aj^) 
± ^sin29K. sin26TCOS<I>Tlni(AoA_L) 

Normalization of the distribution to unity implies |^oP + + I^J-P = 1; and an unobservable phase is removed 
by requiring arg(^o) — 0- This leaves four measurable quantities: the polarization fractions, Fl/F = \Aq\'^ and 
r±/r = |^_Lp, and the phases of the matrix elements, arg(A||) and arg(ylj^). The last two terms of the angular 

distribution have opposite signs for particle and anti-particle decay. The B° and B° decays are flavor tagged by 
the charge of the K meson. The B^ and B^ are not distinguishable by their final state particles, so for B^ decays 
information about the phase of A± is lost. 

The matrix elements are extracted by fitting the observed decay angular distribution. The fit is performed over 
the decay angles and the B candidate mass range covered in Figure 1. To account for detector acceptance and 
selection criteria, Monte Carlo events are generated with detector and trigger simulations and processed by the same 
reconstruction software used on the data. The fit method uses an unbinned log-likelihood with the normalization 
depending on the parameters in the fit [|l^,^ . 

The background is modeled by a sum of polynomial terms of cos Qk* and cos Ot, and sines and cosines of <i>T and 
2<i>x. The background angular distribution is determined from the events on both sides of the B meson mass peak. 

Additional terms are included in the angular distribution to account for the residual probability ('--'6%) of misidenti- 
fying the final state hadrons in the selected B^ J/ip K*^ decays. Events reconstructed under the wrong hypothesis 
will yield incorrect decay angles having a distribution different from, but fully correlated with, the polarized decay 
distribution. The augmented likelihood function corrects for the effect based on the kinematics as determined from 
Monte Carlo. 
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The results from the global fit are illustrated in Figure y. The independent variables are functions of the transversity 
angles: cos Qk*, cos 6t, and $t- The points are the background subtracted projection of the data corrected for 
detector and reconstruction efficiencies. The superimposed curves, derived from the results of the global fit, are in 
good agreement with the data. 

Sources of systematic uncertainty are, in order of importance, the model of the background shape, Monte Carlo B 
generation parameters, trigger simulation, and B^-B^ cross talk. The last two are found to have negligible contribution. 
The uncertainty of the final result is still dominated by the statistics of the event sample [Q. 

The underlying background shape is not known from first principles. To estimate the systematic uncertainty from 
the background modelling, multiple fits with different parametrizations of the background decay angles were studied. 
Different models of the shape are sensitive to different inaccuracies of the background shape function, and the spread 
in the fitted values is used to estimate the sensitivity to possibly unaccounted structure. 

The input parameters to the Monte Carlo generation of B mesons are the h quark mass and generation mass scale, 
the parton distribution functions and fragmentation parameters. The effect on the measurement is determined by 
varying the parameters by their nominal uncertainties and refitting for the decay amplitudes. Systematic uncertainties 
due to the modeling of the trigger selection criteria are studied in a similar fashion. 

The similarity of the decay kinematics for the B^ and B^ decays allows for a possible cross contamination through 
K-TT misidentification. Monte Carlo studies give contamination fractions of 2.1% (4.4%) in the {B^) sample. The 
assigned systematic uncertainty is equal to the difference in the B^ and 5° polarization fraction, multiplied by the 
contamination fraction. 

The final results of this analysis are summarized in Table |, and are illustrated in Figure ||. For the i?" J/ip K*'^ 
the magnitudes of decay amplitudes are in good agreement with the results from CLEO The longitudinal 

polarization fraction dominates: Fl/F = |^oP = 0.59 ± 0.06 ± 0.01 and the parity-odd fraction is small: Fj^/F — 
\Aj_\^ = O.lS^pgg ± 0.06. In contrast to the CLEO result, the CDF measurement for the phase of Ay favors 
the presence of strong final state interactions, in as much as arg(A||) is not close to or tt. However, the results 
from the two experiments are not inconsistent, given their uncertainties. The results for the B^ —> are the 

first available from a full angular analysis. The B^ results also show a dominant longitudinal polarization fraction 
Fl/F = |Aop 0.61 ± 0.14 ± 0.02 and a small parity-odd fraction: F_l/F = = 0.23 ± 0.19 ± 0.04, consistent 

with the B^ results and SU(3) flavor symmetry. 
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Qua.iit ity 


B" .l/rK*" 


B'J J/lu,) 


Ao 


0.77 ±0.04 ±0.01 


0.78 ±0.09 ±0.01 




0.53 ±0.11 ±0.04 


0.41 ±0.23 ±0.05 


arg{A\\) 


2.2 ±0.5 ±0.1 


1.1 ± 1.3 ±0.2 


Al. 


0.36 ±0.16 ±0.08 


0.48 ± 0.20 ± 0.04 


arg{A±) 


-0.6 ±0.5 ±0.1 




Fl/F = 


0.59 ±0.06 ±0.01 


0.61 ±0.14 ±0.02 


Fx/F = 


0.13+|';Jq ±0.06 


0.23 ±0.19 ±0.04 



TABLE I. Fitted decay amplitudes and phases (in radians). The uncertainties are statistical and systematic, respectively. 
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FIG. 1. Invariant mass distributions, after all selection requirements, for (a) ij" ^ J/ipK*^ and (b) J/4'(t> candidates. 



(a) B"^J/vK* 




7C -7t/2 7t/2 It 



FIG. 2. Projections of the results from the full angular fit and the background subtracted data, corrected for acceptance, 
are shown for each of the decay angles, for (a) the J/ip K*'^ mode and (b) the J/ip (f) mode. 



(a) B" ]/\f K* 









II 




FIG. 3. One standard deviation contours, including statistical and systematic uncertainties, for the fitted decay amplitudes 
in the complex plane, for (a) B^' J/ip K*'^ and (b) B^ J/'4>4' decays. 
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